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Abstract 



We study the partition function of both Close-Packed Dimers and the 
Critical Ising Model on a square lattice embedded on a genus two surface. 
Using numerical and analytical methods we show that the determinants of the 
Kasteleyn adjacency matrices have a dependence on the boundary conditions 
that, for large lattice size, can be expressed in terms of genus two theta 
functions. The period matrix characterizing the continuum limit of the lattice 
is computed using a discrete holomorphic structure. These results relate in 
a direct way the lattice combinatorics with conformal field theory, providing 
new insight to the lattice regularization of conformal field theories on higher 
genus Riemann Surfaces. 
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I. INTRODUCTION 



When KaufmanS first evaluated the finite size partition function for the Ising model on a 
toroidal square lattice, she found that the trace of the transfer matrix could be expressed in a 
very compact form as the sum of four terms. In later solutions of the same problem, both in 
the combinatorial approaches!!! or in the closely related Grassmann variable approachS, the 
sum of four terms appears as a natural way of expressing the torus partition function. For 
the Close-Pack Dimer problem, on the same lattice, Kasteleyn^l showed that the partition 
function is given also by the sum of four terms, one of them found to vanish after explicit 
calculation. 

Kasteleyni'ii developed a combinatorial method of solutionEHl that can be applied to 
both the Ising model and the Dimer problem. In this method the four terms are the Pfaffians 
of four adjacency matrices corresponding to different lattice edge orientations 



Z g=1 (T) = - (-Pf (A,) + Pf (A 2 ) + Pf (A 3 ) + Pf (At)) • 



The Pfaffian of an anti-symmetric matrix is the square root of the corresponding deter- 
minant. These determinants can be explicitly evaluated for a square lattice with M rows 
and N columns. For the Ising model at criticality and for the Dimer problem with any dimer 
weights, one of the Pfaffians vanishes and the three remaining behave asymptoticallytiittJ, in 
the M = MN — > oo limit with a fixed ratio M/N, as a common bulk term times a Riemann 
theta functional fc (O, r) of even characteristic. In this limit, the partition function for both 
models can be written as 



e 2 (o\r) 



t](t) 



3 (O|t) 



+ 



4 (O|r) 



exp(/ d A0 (l + O[(logA0 3 /A/]) (1.2) 



with d = 1, 2 for Ising and Dimers respectively, /<j being the free energy of the model. The 
modular parameter in the theta functions given by 



%M cosh 2Kf l /N cosh 2K£ for Ising 



iMz h /Nz v 



for Dimers 



;i.3) 



with K%//3 and K%/ (3 being the vertical and horizontal coupling constants of the Ising model 
at the critical point and z v and z^ the vertical and horizontal Dimer weights. 

Equation ( |1.2[ ) provides a bridge between lattice combinatorics and conformal field the- 
ory. For the Ising model the term in square brackets is the modular invariant partition 
function of the c = \ conformal field theory on a torus, while for Dimers it correspond 
to a c = 1 conformal field theoryS. The theta function dependence of the determinants 
of the different adjacency matrices reproduce the dependence of the determinant of the 
Dirac operator over the different boundary conditions (or spin structures) of the conformal 
field theory&§. This kind of relation between combinatorics and analysis is well know in 
the mathematics literature in the context of the Ray-Singer theoremBEl and it has been 
mentioned in the context of the Ising model on the torus by Nash and O 'Connoi@§. 



In this paper we study the extension of (|1.2| ) to higher genus lattices. For a lattice embed- 
ded in a genus g surface, the partition function of both Ising and Dimers can be expressed as 
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the sum over the Pfaffians of 4 9 adjacency matrices!'!!!. The number of different Kasteleyn 
orientations 4 5 being precisely the number of different fermion boundary conditions in the 
corresponding genus g Riemann surface. We expect the higher genus lattice determinants 
to be related to the functional determinants of the conformal field theory in a way similar 
to the toroidal case. 

While expression (|1.2|) has been known for the torus for over thirty years, the corre- 
sponding higher genus case has not been studied. Simple non-orientable topologies have 
been considered, such as the Mobius stripS'll and the Klein bottle§, and lattices with spe- 
cial symmetriescl including the genus three Klein group lattice L(2,7)II, a lattice with 168 
vertices where the thermodynamic limit cannot be considered. Three dimensional lattices 
have also been studied with the Kasteleyn formalism by embedding the three dimensional 
lattice on a two dimensional lattice of genus of the order of the lattice sizell. While the 
Kasteleyn formalism is well understood for any lattice embedded on an orientablei or a 
non-orientablei surface, explicit computations of the adjacency matrix determinants are not 
known for higher genus lattices. 

We investigate these problems by studying Dimers and the critical Ising model on the 
genus two lattice shown in Fig. |l|. We show numerically that the corresponding sixteen 
Pfaffians of adjacency matrices converge to genus two theta functions and that in the large 
M limit the partition functions of both models can be written as 



Z 9 =2{T C ) 



16 



.1=1 



A exp(f d M) (1-4) 



where d — 1,2 for Ising and Dimers respectively. The 0[z](z|f2) are the sixteen genus 
two theta functions^! with half-integer characteristics, defined in section |TTT| . The factor A 
includes normalization factors and non bulk terms and is not obtainable by our method 
that consists in the numerical evaluation of ratios of determinants of adjacency matrices for 
large lattice sizes. The term in square brackets reproduces the classical winding part of the 
corresponding conformal field theory partition function00. 



While on the torus, the modular parameter (|1.3|) can be understood as a weighted lattice 
aspect ratio, on higher genus lattices the dependence of the period matrix Q of equation ( jl.4| ) 
on the lattice properties is more elusive. We will study this dependence by introducing a 
theory of discrete differentials on the lattice and extract the period matrix from solutions of 
the finite difference Poisson equation on that lattice. The mathematical motivation for this 
procedure is the idea of approximating smooth objects defined on a surface by combinatorial 
objects defined on a triangulation of that surfaceS'ljI'liE. 

The period matrix characterizing a surface is evaluated in two steps, the first step is the 
determination of the space of harmonic differentials on that surface, the second step being 
the decomposition of this space into the holomorphic and anti-holomorphic sub-spaces. 

Harmonic differentials are determined by Hodge decomposition^!! starting from an inner 
product defined on differential forms. On a lattice embedded on a surface the role of p-forms 
is played by the p-cochains, which are linear functionals defined on formal sums of the lattice 
p-elements: vertices, edges and faces. We will call them respectively the lattice functions, 
lattice differentials and lattice volume forms. Eckmannil showed that any choice of inner 
product on the p-cochains gives rise to a combinatorial Hodge theory and to a definition of 
harmonic p-cochains. Later Dodziukil and Dodziuk and Patodi0 proved that for a suitable 
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choice of this inner product on p-cochains, of a triangulation of a surface, the resulting 
combinatorial Hodge theory converges, in the continuum limit, to the surface Hodge theory. 

While our definitions are based on the same general idea, replacing p-forms on a surface 
by p-cochains on a lattice embedded on that surface, we depart from Dodziuk and Patodi 
work on several points. The main difference is that their definition of the inner product on 
p-cochains depends on a mapping of these cochains into forms on the embedding surfaceii 
while our definition is done in term of lattice quantities. Dodziuk and Patodi discrete theory 
is defined referring to the surface properties while we are interested in the inverse procedure: 
to built the continuum theory from the combinatorics of the discrete lattice. 

To determine the period matrix we need the additional concept of discrete holomor- 
phic differentials. These problems have been considered in the study of finite element 
electro dynamicsfEHll where it is well known that the concept of holomorphic differentials 
"defy a straightforward discretization" EH. In the context of the Ising model, discrete holo- 
morphy can be traced back to the finite difference equations on the correlation functions 
found by McCoy, Perk and Wu@ and has been recently discussed by Mercatffl. In Mercat's 
work the difficulty of defining discrete holomorphy for a finite size lattice is patent from the 
fact that the various spaces of differentials have dimensions double of the continuum ana- 
logues. 

In this paper we give a construction procedure for holomorphic differentials, based exclu- 
sively on the direct lattice, which becomes meaningful for large lattice size. To the numerical 
precision we were able to test it, this procedure has a well defined continuum limit and repro- 
duces the period matrices observed in equation (|1.4j ). This is the first time that an argument 
is given to evaluate modular parameters for higher genus lattices. 

The structure of the paper is as follows: 

In section [Djwe review the Kasteleyn formalism and specify the adjacency matrices for 
both the Ising model and the Dimer problem. In section |Hj we establish (|1.4j) by showing 
numerically that ratios of determinants of the adjacency matrices converge to ratios of theta 
functions. In section |Iy| we define discrete harmonic differentials on the lattice and give 
a procedure for numerical evaluation of the corresponding basis. In section ^ we discuss 
discrete holomorphic differentials and the computation of approximations to the period 
matrix. In section [VI| these evaluations of lattice period matrices are compared with the 
period matrices of ( |1.4| ), obtained by direct fit of determinant ratios. Our conclusions are 
presented in section |VII| and the technical details of the Kasteleyn formalism are discussed 
in appendix 0. 



II. DETERMINANTS OF ADJACENCY MATRICES 

Consider the square lattice with the boundary conditions shown in Fig. |l]. This lattice, 
G, is characterized by five integers sizes (Mi, M 2 , K, Ni, N 2 ) and two coupling constants or 
dimer weights (zh,z v ) distinguishing the vertical and the horizontal directions. The boundary 
conditions are such that the lattice can be drawn without superposition of edges only on a 
surface of genus two or higher. The cycles of lattice edges cij,&i, with i = 1,2, represent a 
canonical basis of the first homology group of the embedding surface. The same basis drawn 
over the dual lattice edges will be denoted by a^, 6j (see Fig. 0). 
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We start by considering the Close-Packed Dimer problem on this lattice and for simplicity 
we assume that all integer sizes (Mi, M 2 , K, Nx,N 2 ) are even integers. A close-packed dimer 
configuration in G is a selection of edges such that every vertex is included once and only 



once as a boundary of an edge, see Fig. [T2| for examples. If a weight z h is assigned to the 
horizontal edges and a weight z v to the vertical edges, the dimer partition function is defined 
to be 

G 

z Dimer = £ ^ 
dimer config. 

where and n v are the number of horizontal and vertical edges in a given dimer configu- 
ration and the sum runs over all possible dimer configurations on G. 

While the Close-Packed Dimers problem is an interesting statistical mechanics model 
on its own it also provides a combinatorial approach to the Ising model. There is a well 
known correspondence between polygon configurations of the Ising model high temperature 
expansion and dimer configurations on a decorated lattice!. The Ising model partition 
function on the lattice G can then be expressed as 

G> 

Rising = ( 2 cos h cosh J2 Wv nv w h nh (2.2) 

dimer config. 

with Wi = tanh /3i£j, where K v and K% are the vertical and horizontal coupling constants and 
M is the number of vertices in the lattice G. The sum runs over all the dimer configurations 
of the decorated lattice G' represented in Fig. |. 

Kasteleyn developed a combinatorial formalismiiHi that allows the expression of these 
dimer partition functions as a sum over the Pfaffians of 4 9 adjacency matrices, where g is 
the genus of the simplest surface where that lattice can be drawn without superposition of 
edges. A detailed discussion of the higher genus Kasteleyn formalism is given in appendix 
[A], here we will state only the main results. 

The Kasteleyn adjacency matrices are defined in the following way: label the lattice 
vertices with an integer from 1 to M and choose a lattice edge orientation by assigning to 
each edge a direction represented by an arrow (see Figs. § and |] for two examples). The 
signed adjacency matrix corresponding to this edge orientation is the M x jV matrix Aij 
with entries 

!z if there is an arrow from vertex i to vertex j of weight z 
—z if there is an arrow from vertex j to vertex i of weight z (2.3) 
otherwise 

The Pfaffian of such an JV x M anti-symmetric matrix, with M even, is defined as 

Pf (A) = 2J\T/2^j\f ^2)1 51 €p ^Pip^P3P4 ' ' ' Aw-ipaa (2-4) 

where the sum goes over all the permutations p of the integers from 1 to M and e p = ±1 for 
even and odd permutations respectively. 

From the definition of the adjacency matrix it is clear that each non zero term in the 
Pfaffian expansion equals, in absolute value, a term in the corresponding dimer partition 
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function. The relative sign between different terms depends on the choice of the edge 
orientation. We would like to choose an orientation such that all the terms in the Pfaffian 
have the same relative sign. Then the Pfaffian of the adjacency matrix would equal the 
partition function modulo an overall sign. 

Kasteleyn showed in a beautiful tour de force of combinatorics^ that edge orientations 
with this property exist. These are the edge orientations such that every lattice face was 
an odd number of clockwise oriented edges. In a genus g lattice there will be 4 9 different 
such orientations and the dimer partition function is given as a linear combination of their 
Pfaffiansiii. 

For the genus two lattices G and G', that we are considering in this paper, the relevant 
sixteen edge orientations can be labeled as nr , na 2 , nr ) with n x = 0,1 for x = 

ai,a 2 ,bi,b 2 . The starting orientations v4(0000) are shown in Fig. [| and Fig. |] for the 
lattices G and G' respectively. An orientation with a certain n x = 1 is obtained from the 
corresponding orientation with n x = by introducing a disorder loop along the non-trivial 
cycle x, this is, by reversing the orientation of all the edges crossed by the cycle x. 

To make connection with the theta functions characteristics and allow for more compact 
equations we will also use the alternative notations 



A 



bo 
l d-2 



A; 



(2.5) 



4n~ bi - 2ria 2 - n~ b2 . 



with the integer label given by % = 16 — 8n- c 

We show in appendix |S| that, in terms of these orientations, the Dimer and the Ising 
model partition functions on G are given by 



Pi ' P2 

Pg + Pio 



Ps ~ Pa 

Pll + Pl2 



Pr> 
Pl3 



Pg 

Pl4 



+ P 7 + P S ~ 
+ Pl5 + Pw ) 



(2.6) 



where Pi = Pf Ai and the adjacency matrices are defined on the lattice G for Dimers and 
on the lattice G' for the Ising model. 

For translational invariant lattices, as the torus square lattice, the Pfaffians of the adja- 
cency matrices can be evaluated in a closed form and the theta function dependence can be 
extracted by a careful asymptotic analysis0Hll. For the genus two lattices G and G' such 
an analytic treatment is not possible and we are forced to resort to numerical evaluations of 
the Pfaffians. 



III. RATIOS OF DETERMINANTS AND THETA FUNCTIONS 

Pfaffians of adjacency matrices can be numerically evaluated for large lattice sizes us- 
ing the fact that Pf A = \J det A for an anti-symmetric matrix A. For a given lattice as- 
pect ratio (rrii, m 2, k, n\, n 2 ) we evaluate the ratios of determinants of the adjacency ma- 
trices for a sequence of lattices characterized by the integers sizes (M 1 , M 2 , K, N 1: N 2 ) = 
(mi,m 2 , k,ni,n 2 ) L with increasing L. The Ising model coupling constants were chosen to 
satisfy the criticality condition for a square lattice on the torus 

smh(2f3K h ) sinh(2/3K v ) = 1 (3.1) 
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while no restrictions were imposed on the vertical and horizontal weights of the Dimer model. 
Numerical examples for the critical Ising model and Dimers are shown in tables [j] and ||. 

The fifteen ratios of determinants, obtained is this way, are found to converge smoothly 
with the lattice size TV, see Fig. |], || and ^. The solid line in Fig. |^ and || is a fit with a 
quadratic polynomial on 1/jV to the values obtained for different lattice sizes. 

Unlike the toroidal case, none of the sixteen determinants vanishes at finite size for 
the ( |3.1| ) choice of the Ising coupling constants or for any choice of vertical and horizontal 
dimer weights of the Dimer problem. There are choices of Ising coupling constants with 
small deviations from the criticality condition ( |3.1| ) that will make some of the determinants 
vanish at finite size but not all of them simultaneously. However six of the determinant 
ratios converge to small values for large M and can be associated with theta functions of 
odd characteristic. The convergence of these six ratios is shown in Fig. |7| where the values 
corresponding to the Ising model ratios are squared for comparison. 

The genus two theta functions G are defined by the quickly converging series expansion^ 
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ct 



(z,0) 



0) 



(3.2) 



exp 27r(n + a) T 0(n + a) + 2ni(n + a.) 

n€Z 2 

where a, (3, z and n are 2-vectors half-integers, complex numbers and integers respectively. 
The 2x2 period matrix O is a symmetric complex matrix with positive definite imaginary 
part. 

We found that the extrapolated values of the determinants ratios, in the L 
can be expressed in terms of ratios of theta functions as 

det [A 



oo limit, 



det (A 
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(3.3) 



for the 16 combinations of q = 0, 1 with d = 2 for critical Ising and d — 4 for Dimers. The 
period matrix O being determined from the determinant ratios by a suitable numerical fitting 



procedure. See tables |T| and |V| for examples. Each table displaying the results for three 
different lattices, the first column shows the L —>■ oo extrapolated ratios of determinants 
and the second column shows the theta ratios 



(16-8di -4ci -2d 2 -ca) (^) 







ci/2 c 2 /2 
di/2 d 2 /2 



(o,o) /e ° ; (o,o) 



(3.4) 



for q, dj = 0, 1, with a period matrix obtained by numerical fit. The precision to which the 
two sets of numbers agree is remarkable, typically a precision from 10~ 4 to 10 -6 . 

For the Ising model, way from the criticality condition ( |3.1| ), the sixteen determinants 
are found to converge rapidly to the same value, while for the Dimer problem the theta 
function expression ( |3.3| ) is valid for all values of the dimer weights. 

The period matrices found in ( |3.3|) are purely imaginary for all lattice aspect ratios and 
coupling constants. For both Ising and Dimers we have 

On O12 
O12 O22 



(3.5) 



For such a period matrix the theta functions at zero argument are real. This seems to be a 
property of locally square lattices observed on the toroidal square lattice ( |1.3| ) but not on 
the corresponding triangular latticeEl where the modular parameter is in general complex. 
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IV. HARMONIC DIFFERENTIALS ON THE LATTICE 



In this section we will consider quantities defined on the lattice p-elements: vertices, 
edges and faces. We will call them respectively the lattice functions, lattice differentials and 
lattice volume forms. 

A lattice function / is defined by its value at each lattice vertex and can be represented 
by a M- vector f[n]:n = l,...,Af after some integer labeling of the M vertices of the lattice 
is chosen. 

A lattice differential w is defined by its value on the oriented lattice edges and will be 
represented by a M x 2-matrix io[n|p] : n = 1, . . . , A/"; p = 1,2 where [n|l] stands for the 
edge immediately right of the vertex n and [n|2] for the edge immediately above, referring 
to the lattice drawing of Fig. |l[ We define the horizontal edges to be oriented from left 
to right and vertical edges from bottom to top. The line integral of a lattice differential w 
along a path C of lattice edges is the sum of the values of w on all the edges included in the 
path 

w = ±w[ n \p] (4- 1 ) 

C Hp]ec 

the minus sign applying to edges with opposite orientation to the one of the path. 

A lattice volume form 77 is defined by its value on each lattice face. The lattice in 
Fig. |l] when drawn on a genus two surface has M — 2 faces, two of which are octagons 
the remaining being squares. The volume form 77 will be represented by a (A/"-2)-vector 
r)[q] : q = 1, . . . , Af — 2 where q is an integer labeling the M-2 faces of the lattice. The 
integral of a volume form 77 over a given area A on the lattice is the sum of the values that 
77 takes on all the lattice faces contained in that area 



IL 



v=E v[q] (4-2) 

A [<AeA 

It is convenient to relate the labeling of vertices with the labeling of faces. For that 
purpose we introduce the following notation: if q labels a squared face then q\ stands for 
the label of its lower left vertex; if q labels an octagonal face then q\ and q\ stand for the 
labels of the two vertices that can be seen as the octagons lower left vertices in Fig. [T|. 
Reciprocally if n labels a lattice vertex then n is the label of the lattice face of which n can 
be seen as the lower left vertex. We then have the relations 

ii = q (4.3) 

n e {n 1: n 2 } (4.4) 

An inner product on lattice functions, lattice differentials and lattice volume forms can 
be defined in the following way 

(/,/0=E/W/ / W ( 4 - 5 ) 

n 

[W, W ) = (h/v) w[n\l]w'[n\l] + {h/v)~ l w[n\2]w'[n\2}) 

n 
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where h/v is a positive parameter providing different weighting of the horizontal and vertical 
components of a differential. 

The lattice exterior derivative d is a linear operator that takes lattice functions into 
lattice differentials and lattice differentials into lattice volume forms. For our labeling and 
choice of orientations d is given by 

(d/)[n|l] = /[right(n)]-/[n] (4.6) 
(df) [n\2] = /[up(n)] - f[n] 
(dw) [q] = J2 (Mftll] + Mright(gi)|2] - ™[up(^)|l] - ^|2]) 

i 

with i=l,2 for octagonal faces and i—1 for squared faces. The functions 
right(n), left(n), up(n) and down(n) give the label of the vertex immediately right, left, above 
and below of vertex n, referring to the lattice drawing of Fig. |]. 

The exterior derivative defined in this way satisfies a discrete version of Stokes theorem. 
Let C (n, n') be a path of lattice edges from vertex n to vertex n' and A an area on the 
lattice, then we have that 

df = f[n']-f[n] (4.7) 

C(n,n') 

dw = w 

JdA 

where dA is the path along the boundary of the area A with an anticlockwise orientation. 
The co-derivative 5 is defined as the adjoint operator of the exterior derivative 

(w,df) = (5w,f) (4.8) 
(rj, dw) = (5r), w) 

and can be expressed as an operator that takes lattice differentials into lattice functions and 
lattice volume forms into lattice differentials 

(5w) [n] = - (w[left(n)|l] - w[n\l}) + ~ (tu[down(n)|2] - w[n\2]) (4.9) 

v / — 
(5r)) [n\l] = — [T][n\ — r][dovjn(n)} 

(5 7]) [n\2] = ^ (?7[left(n)] - rj[n 



(4.10) 



where left(ra) stands for the face given by the tilde of the vertex left(n) and similar for 
down(n). 

Both the exterior derivative and the co-derivative satisfy dd = 55 = 0. Explicitly in 
terms of the formulas above we have that 

(ddf) [q] = (/[up(rightfo))] - /[right(up(ft))]) = 0, (4.11) 

i 

(55 rj) [n] = r;[left(down(n))] - r/[dovvn(left(n))] = 0. (4.12) 
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There are four vertices on the octagonal faces where right(up(n)) 7^ up(right(n)), at these 
faces dd vanishes because of cancellation between the qi and the q% terms. There are also 
four vertices where left(down(n)) 7^ down(left(n)) but left(down(n)) = down(left(n)) is true 
for all vertices. 

In exact analogy with the continuum definitions, a lattice differential w is said to be 
closed if dw = 0, exact if w = df, co-closed if 5w = and co-exact if w — 5r). We are 
mainly interested in the harmonic differentials, a lattice differential is said to be harmonic 
if it is both closed and co-closed 



In general, on a genus g lattice, the lattice harmonic differentials form a vector space 
of dimension 2g. For the genus two lattice G this can be seen directly from the system 
of equations ( }4.13j ) . It is a system of 2j\f — 2 equations on 2Af unknowns and there are 
at least two independent solutions. Two additional solutions are provided by constant 
differentials with independent vertical and horizontal components and the system has at 
least four independent solutions. 

It is remarkable that, on the lattice, the dimension of the space of harmonic differentials 
is a linear algebra problem determined by the Euler characteristic. Consider a genus g 
lattice with the property that its edges can be distinguished into two classes: vertical and 
horizontal. For such a lattice the dimension of the space of harmonic differentials, that is, 
the number of independent solutions of the system (|4.13|) on that lattice, is at least 



where \ is the Euler characteristic of the lattice. 

It still remains to be proven that there are not more than 2g linear independent harmonic 
lattice differentials. Note that 2g is the dimension of the lattice first homology group. We will 
prove, in a way similar to the continuum, that a harmonic lattice differential is completely 
determined by its integrals, or periods, around the lattice non-trivial loops and therefore 
there can only be 2g independent harmonic differentials. 

Harmonic differentials are by definition closed and its integrals depend only on the ho- 
mology class of the path considered. If two differentials w and w' have exactly the same 
periods along the 2g classes of loops then their difference w — w' has zero integral along any 
closed path on the lattice and the function 



is well defined for any lattice path going from a fixed vertex o to the vertex n. This function 
is a harmonic function, in the sense that it satisfies 




(4.13) 



#(edges) - #(vertices) - #(faces) + 2 = ~x + 2 = 2g 



(4.14) 




(4.15) 



A/[n] = (5d f) [n] =S(w- w') [n] = 



(4.16) 



with the Laplacian operator given explicitly by 



(A f)[n] = 2(h/v + v/h)f[n] 

-V«(/[nght(n)]+/[left(n)]) 
-v/h (/[up(n)] + /[down(n)]) 



(4.17) 
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It is easily seen that a function harmonic at n cannot have a local extremum at n and that 
the only functions harmonic everywhere on the lattice are constants. From this it follows 
that the function defined in Q4.15Q vanishes at all lattice vertices. In this way, we proved that 
if two harmonic differentials w and w' have the same periods then they are equal, w = w'. 
The space of harmonic differentials has therefore the same dimension as the first homology 
group, 2g. 

For a given lattice a basis of the space of harmonic differentials can be evaluated by 
direct solution of the linear system of equations ( }4.13| ). Such a solution can always be 



found, at least numerically, although the evaluation of the kernel of a large linear system is 
computationally demanding. For practical purposes there is a much more efficient method 
based on a discrete version of the Hodge decomposition theorem. 

The discrete Hodge decomposition theorem states that the space of lattice differentials 
has an orthogonal decomposition in terms of exact, co-exact and harmonic differentials and 
that any lattice differential w can be written in an unique way as 

w = df + h + 5r] (4.18) 

where h is a harmonic differential. The orthogonality of the different kinds of differentials 
follows directly from ( |4.13| ) and the property dd = 55 = 0. 



This result allows the determination of harmonic differentials by orthogonal projections. 
Our objective is to obtain a basis {A k , B k ■ k — 1, 2} of the space of harmonic differentials 
satisfying the normalization conditions 



A k = 5 kj , A k = (4.19) 
B k = 0, B k = 5 kj 

with the a,j,bj being any choice of closed paths on the lattice representing a basis of the 
first homology group. We can proceed in the following way: start with closed but not 
harmonic differentials A k ,B k with the periods required by ( |4 . 1 9| ) . A possible choice for 



these differentials is shown in Fig. |8| where we take A k (B k ) to be zero on all edges except 
the ones crossed by the dual lattice cycles b k (respectively d k ). 

Since closed differentials are orthogonal to co-exact differentials it follows from (|4.18| ) 



that closed differentials can be written as the sum of a harmonic differential with the same 
periods and an exact differential 



A k = A k + dfC (4.20) 
B k = B k + df^ 



applying 5 on the these equations and using fl4.13|) and ( |4.16|) we obtain 



Af£ = 6A k (4.21) 



^fk=SB, 



The Laplacian operator on functions (|4.17|) is a well defined M x jV matrix of rank 
(A/ - — 1) acting on the functions A/"-vector and ( |4.21| ) can be solved numerically for reasonably 
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large M by fixing the value of the functions at a lattice vertex. The solutions can then be 
differentiated and subtracted to the A k ,Bk to obtain a normalized basis of the space of 
harmonic differentials {A^,Bk : k = 1,2}. This method is numerically more efficient than 
the direct solution of the linear system ([4.13|) . 



V. HOLOMORPHIC DIFFERENTIALS AND THE PERIOD MATRIX 

To evaluate a discrete approximation to the period matrix we need to introduce the 
concept of lattice holomorphic differentials and decompose the space of harmonic differ- 
entials {Afc, -Bfc : k = 1,2} obtained in the previous section into the holomorphic and 
anti-holomorphic sub-spaces. 

In the continuum theory this is accomplished by the projection operators 

P=(l + i*)/2 and P=(l-i*)/2 (5.1) 



where * is the (continuum) Hodge operatorclH This operator is an endomorphism on the 
space of harmonic differentials and satisfies * 2 = — 1. Its action on differentials is given by 

* (fx(x, y) dx + f y (x, y)dy) = -f v (x, y) dx + f x (x, y)dy (5.2) 

and a differential is said to be holomorphic if it is of the form Pw with w harmonic. 

We will proceed in a similar way on the discrete theory and define lattice holomorphic 
differentials by a projection with a discrete Hodge operator *. It proves very difficult to 
define a discrete Hodge operator acting on lattice differentials, on the finite size lattice, with 
the same properties of the continuum one. We propose the following definition 

(*w)[n|l] = -w[down(n)|2] (h/v)' 1 (5.3) 
(*w)[n|2] = w[left(n)|l] {h/v). 

The action of the discrete Hodge star on the lattice differentials is shown in Fig. |9|. In 
this figure an arrow connecting two edges represents that the value of *w at the second edge 
is obtained from the value of differential w at the first edge. Besides the definition above 
there are three related definitions of the Hodge star, corresponding to it/2 rotations, that 
are also shown in Fig. ^. Due to the symmetries of the lattice, the four different definitions 
produce the same numerical results for the period matrix. 

The principal merit of these definitions is the way they relate the exterior derivative d 
with the co-derivative 5. From the definition above and ( |4.6|) and ( [4.9D it follows that 

{d*w)[q] = - y £6w[q i \ (5.4) 

i 

j a(n) if n e {ni,n 2 ,n 3 ,n 4 } 

[d*w)\n\ = { . r . — /, r / w-i i • (5-5) 

v n 1 \ dw[down(left(ra))] otherwise v ' 

where in the first equation the sum goes over i=l,2 for octagonal faces and is i=l for squared 
faces. The are the four lattice vertices where left(down(rafc)) ^ down(left(n&)). Let n\ 
and n 2 be on the same octagon and the two remaining points on the other octagon then the 
a{n k ) satisfy 
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a(ni) + a{n 2 ) = diy[down(left(ni))] (5.6) 
a(n 3 ) + a(n 4 ) = d u?[down(left(n 3 ))]. 



The appearance of the a terms in ( |5.5|) prevents the discrete Hodge operator from being 
an endomorphism on the harmonic lattice differentials. These terms originate from the fact 
that the lattice has more vertices than faces: the requirement that a differential w is closed 
is not sufficient to ensure that the corresponding -kw is co-closed. For the dual lattice, where 
the number of faces is larger than the number of vertices, the opposite happens and the a 
terms appear in Q5.4D instead of (|5.5|) . 

Unlike the continuum Hodge star, our lattice definition does not satisfy * 2 = — 1 but 
rather 

* 2 w[n\l] = -tu[left(down(n))|l] (5.7) 
* 2 w[n\2] = -tu[down(left(n))|2] 

there is a displacement of the differential and mixing at the octagons. 

Notice that for a toroidal square lattice, where all faces are squares and 
left(down(n)) = down(left(n)) for all vertices n, there are no a(n) terms on the equivalent of 
equations ( |5.4| ) and (|5.5|). The lattice Hodge star ( |5.3j ) on a torus is an endomorphism on 
the space of harmonic differentials. 

Lets proceed to evaluate the torus modular parameter r and compare with the know 
result ( |1.3| ). Consider a toroidal square lattice with M rows and N columns. Let {a, b} be a 
basis of the first homology group with a being a horizontal loop and b a vertical loop. This 
lattice has two independent normalized harmonic differentials A and B with components 

A[n][l] = l/N, A[n)[2) = , , 

B[n)[l] = 0, B[n)[2) = l/M { } 

for all vertices n. The only holomorphic differential on the lattice is V = A + i * A, in 
components, 

r[n][l] = ^, r[n][2]=i^i (5.9) 

for all vertices n. The modular parameter is the integral of this differential along a loop of 
the homology class of b and we obtain 

t= f(A + i*A)=i^ (5.10) 
Jb v IS 

that reproduces the well known result ( |1.3j ) if we take for the h/v ratio 

h/ — f cosh 2^/ cosh 2K% for the critical Ising Model , ^ 

1 Zh/z v for the Close-Packed Dimers 

For the Ising model this ratio corresponds to the ratio of the diverging correlation lengths 
along the two lattice directions^!, the relevant geometrical quantity at criticality. We will use 
the same expression for the genus two case, however in this case the lattice Hodge operator 
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( [5.3| )is not an exact endomorphism on the harmonic lattice differentials, due to the a terms 
that appear on the four n k vertices. 

We will take ( pT3| ) as the definition of our lattice Hodge operator and proceed to evaluate 
a basis for the space of differentials of the form {T^ = (1 + i-k)H k : k = 1 . . . , g} where the 
Hk are harmonic differentials, satisfying the normalization conditions 

/ {H l + i*H l ) = 5 u (5.12) 

[ {H l + i*H l ) = Q kl (5.13) 



Equation d512|) is a set of 2g real constraints on H k from which we can determine the 



coefficients of H k in the basis {A k} B k } and equation ( {5.13Q gives our approximation to the 
lattice period matrix. 

Since {T k = (1 + i*)H k : k = 1 . . . ,g}, seen as a complex vector space, is not invariant 
under P = (1 + i~k)/2 we are not allowed to call these lattice differentials holomorphic in a 
strict sense. However a detailed numerical study, discussed on the next section, will show 
that the period matrix obtained by this procedure reproduces to excellent precision the 
period matrices of equation ([ 



VI. PERIOD MATRICES AND DETERMINANT RATIOS 

Following the procedure described in the two previous sections we can evaluate nu- 
merically the period matrix for lattices with different aspect ratios (m l5 m2, k, n 2 ) and 
coupling constants. Table |V| shows the period matrix for a particular lattice aspect ratio 
(mi, m 2 , k, ni, n 2 ) and various lattice sizes (Mi, M 2 , K, N±, N 2 ) = (mi, m 2 , k, rii, n 2 ) L, with 
increasing L. The period matrix elements converge in a smooth way with the number of 
lattice vertices TV and the resulting matrix is a purely imaginary, symmetric and positive 
definite matrix for all lattice sizes. 

The three different set of values shown in table [V] correspond to three different choices 
of the first homology group basis. The first set of values, labeled as A, corresponds to the 
choice of homology basis Oj, hi shown in Fig. |l|; the set labeled as B corresponds to a choice 
where the ai and hi loops are interchanged from Fig. |T] and the set C corresponds to the 
choice where the top and bottom of Fig. [I] are interchanged. The different period matrices 
obtained for each choice of the homology basis are related by a modular transformation. 
In table |VT| the theta function ratios corresponding to each choice are shown to be related 
by permutations. While the A and C sets of values agree at finite size, the period matrix 
corresponding to B gives numerical values different from the other two. This fact is related 
with the a terms in equations ( |5.4|) and (|5.5| ) that keep the lattice Hodge operator from 



being an endomorphism on lattice harmonic differentials. 

While at finite size, modular invariance is slightly broken this does not seem to be the 
case in the continuum limit. In table |VT|we see that the theta function ratios extrapolated 
to the M — > 00 limit, agree within a 10~ 2 error. Fig. [| and|^ show the theta function ratios 
corresponding to the two different modular choices A and B, plotted as a function of number 
of lattice vertices A/", together with the determinant ratios. In these figures we observe that 
theta ratios for the period matrices A and B seem to converge to a common value and 
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that, for large enough lattice size, the values of the nine non-vanishing determinant ratios 
are bounded by two sets of the theta ratios. Qualitatively the pattern observed in these 
plots leads us to believe that the procedure described on the previous section to evaluate 
holomorphic differentials and the period matrix has a well defined continuum limit. 

A quantitative analysis is given in tables [TTI| and [IV] where the period matrices evaluated 
by the procedure above (using the choice A for the homology basis) are compared with the 
period matrices of equation ( ^.3[ ), obtained from the ratios of adjacency matrices determi- 
nants by numerical fit. There is a remarkable agreement between the two sets of values. 
The theta function ratios corresponding to the two period matrices agree with a precision 
of 10~ 3 or better. 

The small numerical differences are associated with the difficulty in extrapolating the 
finite size values of the period matrix to the continuum limit. The extrapolated values given 
in tables |TTT| , fV[ |V| and |VI] are obtained by fitting the size dependence of the period matrix 



entries with a quadratic polynomial in M . 

As in the toroidal case, the vertical and horizontal coupling constants provide different 
weighting to the vertical and horizontal directions. The period matrix for a lattice charac- 
terized by the aspect ratio (mi, m.2, k, ni, n 2 ) and parameter h/v is equivalent to the period 
matrix of a lattice characterized by the aspect ratio (ami, arri2, (3k, {3ni, fln-i) and a re-scaled 
parameter (3h/av. This property is exemplified in the first two sets of values in tables [Til 
and |Vj for both the critical Ising model and Dimers. The parameter h/v can therefore be 
absorbed into a re-scaling of the overall vertical/horizontal aspect ratio of the lattice. 



VII. CONCLUSIONS 

In this paper we used the Kasteleyn formalism to study Dimers and the critical Ising 
models on a genus two lattice. It is the first time that such a study has been done in 
higher genus lattices. We found that the determinants of the Kasteleyn adjacency matrices 
converge ( j3.3|) to a common term times theta functions of half-integer characteristic. This 
result generalizes a thirtv year old observation on the asymptotics of Pfaffians of adjacency 
matrices on the torusi2rEJ and gives a new meaning to this convergence. 

For the critical Ising model, the dependence of the determinants of adjacency matrices 
on the Kasteleyn orientations is exactly the same as the dependence of the determinant of 
the Dirac operator, of the corresponding conformal field theory, on the spin structures of the 
Riemann surface; both are given in terms of theta functions of half integer characteristics. 
There is therefore an one to one correspondence between the 4 9 Kasteleyn clockwise odd 
orientations of the lattice and the 4 9 spin structures of the continuum limit Riemann surface. 

These observations elevate the Kasteleyn formalism from a combinatorial trick to ob- 
tain the partition function to a discrete formulation of some of the model conformal field 
theory properties. The relation between the lattice determinants and the zeta regulated 
determinants of the conformal field theory is reminiscent of the Ray-Singer theorem@llI 
that relates determinants of the analytical Laplacians on a surface with the determinants of 
combinatorial Laplacians on a triangulation of that surface. 

In the context of conformal field theory the period matrix is usually seen as a free 
parameter and, except on the torus, is not related to the thermodynamic limit of a lattice 
or a set of particular lattices. We have shown that the period matrix, characterizing the 
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continuum limit of a lattice, can be understood in terms of a lattice discrete holomorphic 
structure. 

An important difference between the genus two case ( |1.4| ) and the torus ( |1.2| ) is that the 
determinants corresponding to the odd characteristic theta functions do not vanish at finite 
size and only converge to zero on the M — > oo limit. It interesting to notice that on the 
toroidal square lattice a completely satisfactory definition of discrete holomorphy at finite 
size can be given, a fact that is probably related with the finite size vanishing of the odd 
characteristic determinant. 

The results we obtained for the genus two case can be readily generalized for arbitrary 
genus. Higher genus lattices can be constructed by pasting two or more of the lattices of 
Fig. m, as shown in Fig. |10[ The construction of the adjacency matrices and period matrices 
can also be done in a similar way to the one described in this paper for genus two. 

There are many ways of choosing the boundary conditions of a square lattice to obtain 
a higher genus lattice. An important requirement to reproduce the results we obtained is 
that the boundary conditions do not destroy the distinction between vertical and horizontal 
edges. With the boundary conditions suggested in Fig. [U], the genus g lattice has a total 
of 3g — 1 integer sizes, respectively (Mi, . . . , M g , K\, . . . , K g -x, N\, . . . , N g ). These account 
for 3g — 2 independent aspect ratios, the h/v parameter been absorbed into a combination 
of these ratios. These 3g — 2 rational parameters cannot cover, or even approximate, all 
the complex structures on a genus g Riemann surface, that are parameterized by 3g — 3 
complex parameters. This may be possible by generalizing the genus two lattice of Fig. [I] 
from locally square to locally triangular and allowing for different couplings constants in 
each sub-lattice. 
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APPENDIX A: KASTELEYN COMBINATORICS 



In this appendix we derive ( |2.6j ), the Pfaffian expansion of the partition function. We 
will show that from the Pfaffians of the sixteen Kasteleyn adjacency matrices Pf (Aj) we can 
obtain the partition functions of Ising and Dimers for sixteen different choices of boundary 
conditions on the lattice of Fig. |I[ Let Z(0000) denote the partition function of Dimers or 
Ising on the lattice of Fig. [I], then 

Z t = Z(n dl , n~ bv n~ a2 , ngj (Al) 

with i = 16 — 8na x — 4n^ — 2ria 2 — nr and n x = 0,1 for x = n^^ni^n^^n^ represent the 
partition functions of the same model with coupling constants or dimer weights multiplied 
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by —1 along a choice of the n x = 1 representatives of the first homology group. This is a 
labeling similar to the one introduced for adjacency matrices in ( j2.5|) . 



We will show that there is a matrix 6y such that 

49 49 
Z l = J2 hj Pf (Aj) and Pf (Ai) = £ K Zj (A2) 
3=1 i=i 

The property that b 2 = 1 was first noticed for the torus by Fradkin and Shteingradt@. 
This places the partition functions and the Pfaffians of the adjacency matrices on equal 
footing. For both models the partition function Z(n ai , ng , n a2 , ng 2 ) is independent of the 
choice of representatives of the cycles a&, In the Ising model it is well known that a 
deformation of a disorder loop does not alter the partition function since a disorder loop 
that crosses a vertex corresponds to an interchange of the up and down spin at that vertex. 
For the Dimer problem the invariance of the partition function Z(na u ut, ria 2 , ) follows 
from ( |A2| ) and from the invariance of the determinants of adjacency matrices under disorder 
loop deformation. 

To keep the paper self contained we start by giving a short overview of the Kasteleyn 
formalism in higher genus lattices following a combination of references! - !. By a genus g 
lattice we mean a graph, not necessarily regular, that can be drawn without superposition 
of edges only in a surface of genus g or higher. The Ising model decorated lattice, Fig. |3|, 
has crossing edges and cannot be considered a genus two lattice. To apply the results of this 
appendix to the Ising model it is preferable to consider the Fisher six vertices decorationil 



instead of the Kasteleyn four vertices decoration, see Fig. ^TJ. The Fisher decorated version 
of our lattice can be embedded on a genus two surface without superposition of edges. 
Determinants of adjacency matrices of the two lattices are equali. We will use the Fisher 
lattice on the formal discussion and the Kasteleyn lattice for actual numerical calculations. 
Recall the definition of the Pfaffian of anti-symmetric matrix, equation (|2.4|) 



Pf (A) 2 N / 2 (N/2)\ ^ C ( Pl p ' 2 " PN ) ^ P1P2 ^ P3P4 ' ' ' ^pn-ivn (A3) 

where p is a permutation of the integers from 1 to N and e p — ±1 for even and odd 
permutations respectively . The normalization factor 2 N ^ 2 (N/2)\ accounts for the fact that 
permutations differing by the order of the indices in a A PiPi+1 factor or by the interchange 
of two such factors give equal contributions to the Pfaffian. We will call two permutations 
related in such a way equivalent permutations. 

A non zero Pfaffian term A PlP2 A P3P4 ■ ■ ■ A Pn _ iPn corresponds to a choice of edges C p = 
{(P11P2), • • • , (pn-i,Pn)} where no edges share a common vertex and all vertices are in- 
cluded; (2, j) representing the edge between vertices i and j. This choice of edges is a dimer 
configuration with weight equal, in absolute value, to the Pfaffian term. Dimer configura- 
tions are therefore in one to one correspondence with sets of equivalent permutations with 
non zero Pfaffian term. 

To express the dimer partition function in terms of Pfaffians of adjacency matrices we 
need to choose edge orientations such that all the terms in the Pfaffian have the same 
relative sign or a linear combinations of Pfaffians with this property. The relative sign 
between Pfaffian terms corresponding to different permutations p and p' can studied by 
considering their product 
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e 



( pi P2 ■ ■ Pn) 6 ( P \ k ■ ■ PW A P1P2 ' ' ' A PN-1PN A p[p> 2 ■ ■ ■ A 'p' N ^p' N ( A4 ) 



where A and A' represent the same adjacency matrix, the prime being introduced for conve- 
nience. This product has a simple graphical interpretation: it corresponds to the superposi- 
tion diagram of the two dimer configurations C p and Cy. The superposition diagram is the 
set of double edges and even length cycles obtained by drawing both dimer configurations 
on the lattice, see Fig. O for an example. The even length cycles have alternating edges 



belonging to each one of the two dimer configurations and are called transition cycles. 
The product (S3) can always be rewritten in a form that resembles the superposition 



diagram! 

J *1 *2 ■■ i q iq+l iq+2 \( A. a* A. . A 1 ...A' \( A- A 1 \(...\... ( AK\ 

V *2 13 •■■ h i q +2 iq+l )\^iii2^i 2 i s ^W^i±i5 ^i p i 1 )\^-iq+iiq+2^i q+2 iq+ 1 )\ I \^°) 

where we used the fact that the sign of a product of permutations is the product of the sign 
of each permutation and we decomposed the overall permutation into cyclic permutations of 
even length. We also chose equivalent permutations to p and p' and rearranged the matrix 
elements in such a way that the terms corresponding to the double edges and the even 
length cycles of the superposition diagram are singled out. The equation above is a possible 
example where one cycle of 2p length and a double edge are shown. 

From this formula we can read out the relative sign of two Pfaffian terms. The sign of 
the overall permutation is (— l) s where s is the total number of even length cycles, including 
double edges, in the superposition diagram. From the double edges A^A'^ we get an ad- 
ditional minus sign since A is anti-symmetric, while each transition cycle contributes with 
a sign given by (— l) p where p is the number of edges in that cycle oriented, say, in the 
clockwise direction. 

Following the notation of reference! let C AC denote the set of transition cycles resulting 
from the overlap of the dimer configurations C and C with the double edges removed. Then 
the above discussion motivates the following result due to KasteleynQ 

Theorem 1 Let sign(C) be the sign of the terms in the Pfaffian corresponding to a given 
dimer configuration C. Then for any two dimer configurations C and C we have that 

sign{C)sign{C) = (A6) 

7 

where the product runs over all transition cycles 7 in C AC and ^(7) is the number of 
clockwise oriented edges in the cycle 7. 

We will call A the overlap operator. If dimer configurations are seen as sets of edges 
then 

C 1 AC 2 — (Ci U C 2 )\(d n C 2 ) 

is the symmetric difference of the two sets. Since superpositions of configurations will play 
a major role in the following discussion it is worth to collect some properties of the overlap 
operator A. The following properties follow directly from the definition, the Cj being dimer 
configurations or more generally sets of edges 
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d A C 2 = d A d (A7) 
d A d = 
d A = d 
d A (d A d) = (d A d) A d 

We can use this overlap operator to classify the dimer configurations, and the corre- 
sponding Pfaffian terms, into equivalence classes or types. We start by choosing a standard 
dimer configuration Cq, the standard configuration for our lattice being shown in Fig. [T2. 



Let Oj, bi with % = 1, . . . , g be a canonical basis of the first homology group on the lattice 
then we say that C is of type T(n& 1 , n ai , . . . , nf, g ,n ag ) if there are n x topologically non-trivial 
cycles of kind x = ai, bi present in the overlap Cq A C. An example for our genus two lattice 
is show in Fig. where a dimer configuration C is seen to be of type T(0102). 

We are then classifying C by the homology class of the overlap Cq A C, two dimer 
configurations C and C being of the same type if (Cq A C) = (Cq A C) A X where X is a 
boundary, that is, a set of closed cycles with trivial homology. Since the overlap operator A 
eliminates double edges it can only distinguish homology classes modulo two. For instance in 
Fig. [T3]we show how a T(2001) superposition can be obtained from a T(0001) superposition 
by overlapping a homologically trivial boundary X. We will therefore take n x = e,o where 
were e stands for even and o stands for odd. 

The relevance of these equivalence classes is that the relative sign of a configuration 
C to the standard configuration Cq for a clockwise odd edge orientation depends only on 
the equivalence class of C. Clockwise odd edge orientations were defined in section [Tl] to be 
orientations such that every lattice face has an odd number of clockwise oriented edges along 
its boundary (the clockwise orientation is a convention, anti-clockwise odd edge orientations 
would work equally well). These edge orientations have the following property, found by 
Kasteleyni 

Theorem 2 For an edge orientation such that all lattice faces have an odd number of clock- 
wise oriented edges we have that sign(C)sign(C) = 1 for any pair of dimer configurations 
C and C such that C AC has a trivial homology. 

A complete proof of this theorem can be found in references^ but the general idea is very 
simple: if C A C has a trivial homology its transition cycles can be built by the successively 
overlap of elementary lattice faces with the A operator. It can then be shown that any 
such transition cycle will have clockwise odd parity and contribute with a minus sign in the 
Pfaffian term that cancels the minus sign of the corresponding cyclic permutation. 

We can now state the main result of the higher genus Kasteleyn formalism, whose proof 
is so simple that we present it here 

Theorem 3 The relative sign of a dimer configuration C to the standard configuration Cq 
in a clockwise odd orientation depends only on the homology class modulo two of the set of 
transition cycles Cq AC . 

Proof: Let C and C be two dimer configurations such that Cq AC and Cq A C belong 
to the same homology class. Then by definition we must have (Cq AC) = (Cq A C) A X 
were X is a set of homologically trivial cycles. Using (|A7|) we find that X = C A C and it 
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follows from theorem 2 that, in a clockwise odd orientation, we have sign(C)sign(C > ) = 1. □ 



In a genus g lattice there are A 9 classes of dimer configurations, the number of elements 
of the homology group Hi(G, Z 2 ), that in general will have different signs in a given clock- 
wise odd orientation. Conversely there are also A 9 inequivalent Kasteleyn orientations, the 
number of elements of the cohomology group Z 2 ) and it is possible to find a linear 

combination of the corresponding 4 9 Pfafnans of adjacency matrices such that all dimer 
configurations have the same overall sign. 

More precisely: two edge orientations are said to be equivalent if they are related by a 
sequence of local transformations in which we reverse the edge orientation of all the edges 
coming to a given vertex. It is clear that each local transformation is equivalent to the 
multiplication of the i-th column and row of A by —1 and will not affect the parity ^(7) of 
the transition cycles 7 or the absolute value of det A. We then have the following result 

Theorem 4 In a genus g lattice there are 4 9 inequivalent Kasteleyn edge orientations. 

For our genus two lattice it is easy to see that there are only sixteen inequivalent Kaste- 
leyn orientations. Consider the lattice as drawn in Fig. [14]. We can choose the edge orienta- 
tion on all the edges of a spanning tree of the lattice (heavy line) using the freedom given by 
the local equivalence transformations. Then the condition that the elementary faces must 
be clockwise odd will fix the orientation on most edges (the dashed line edges) but there will 
be 4 edges (rectangular edges in the figure) in which the orientation is undetermined. The 
sixteen different choices of orientations in these edges will correspond to sixteen inequivalent 
Kasteleyn orientations. Using a similar reasoning the reader can convince her/himself that 
the theorem is true for any lattice. 

In section 2 we constructed the sixteen edge orientations for our genus two lattice by 
starting from an initial orientation and introducing disorder loops along the 5j, hi cycles. 
This procedures makes it clear the connection between the edge orientations and the the 
cohomology group i? 1 (G, Z 2 ). 

In order to obtain the linear combination of Pfafnans that yields the partition func- 
tion we still need to determine the relative sign that the different dimer configuration 
types, T(nb 1 n ai n; )2 n a2 ) with n x = e,o, take in the various clockwise odd configurations, 
A{m- a jn- h m- a2 mi^) with m~ x = 0, 1. 

Dimer configurations of type T(eeee) will have the same sign on all clockwise odd orien- 
tations since its overlap with the standard configuration Cq is topologically trivial. The sign 
of a general configuration of type T(nb 1 n ai rib 2 n a2 ) under the edge orientation ^4(0000) can be 
evaluated by inspection on Fig. ^| and f| after choosing a representative overlap diagram of 
each type, remembering that we are considering all the sizes (Mi, M 2 , K, N±, N 2 ) to be even 
integers. This is the last column of table [VI I| . Once these signs are know the remaining signs 
on the table |V11| can be evaluated in a simple way: a relative sign of a configuration changes 



with the introduction of a disorder loop if that configuration has a non-trivial transition 
cycle intersecting that disorder loop. 

Once all the relative signs are known it is simple to find the linear combination of Pfafnans 
that will give the same overall sign to all dimer configurations and we obtain equation ( |2.6| ) 
and the last row of the matrix 6y in equation ( |A2|) . Equation (|2.6|) can be rewritten using 
the notation of (12.51) as 



20 



Z(0000) = ^ £ S^^^Pf^K^^^)) (A8) 

4 n x =0,l 



with some choice of an = ±1. The Pfaffian expansion of the remaining partition functions Zi 
can be obtained form this equation by successive inversions — > 1, 1 — > on the same entry 
n x on both the Z() and the AQ terms and noticing that this correspond to a permutation 
of the aj. In this way we can construct the full matrix bij that for our lattice, and labeling 



of the orientations, is given in table VIII 
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FIG. 1. The genus two lattice with the boundary identifications given by the integers from 1 
to 22. The a^, hi cycles form a basis of the first homology group and should be seen as drawn over 
the lattice edges. All lattice faces are squared except for two octagons whose edges are marked 
with the doted line. 
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FIG. 2. The ^4(0000) clockwise odd orientation of the original lattice. The a>i,bi cycles form a 
basis of the first homology group drawn over the dual lattice edges. The remaining clockwise odd 
orientations can be obtained from this one by inverting the orientations of the edges crossed by a 
given choice of the aj, bi cycles. 
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FIG. 4. The ^4(0000) clockwise odd orientation of the Ising decorated lattice. The remaining 
clockwise odd orientations can be obtained from this one by inverting the orientations of the edges 
crossed by a given choice of the Sj, 6; cycles. 
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FIG. 5. The continuum limit for DIMERS: the nine non- vanishing ratios of determi- 
nants, defined as Ri = det A4/ det Aiq, and the corresponding ratios of theta functions are 
shown as a function of the number of lattice vertices M, for the lattice with aspect ratio 
(mi,m2,k,ni,n2)=(4,2,2,6,2) and 2^=0.780 z„=0.560. The solid curve is the plot of the deter- 
minant ratios and the dashed curves are the plots of ratios of theta functions, with period matrices 
evaluated for the A and B choices of the homology basis (see section VI for discussion) . The curves 
are fits with a polynomial on 1/jV to the values obtained for different lattice sizes. 
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FIG. 6. The continuum limit for ISING: the nine non- vanishing ratios of determinants and 
ratios of theta functions are shown, as a function of the number of lattice vertices jV, for the lattice 
with aspect ratio (mi,m2,k,ni,n2)=(2,2,2,2,2) and u;/j=0.537 u>„=0.301. Some of the ratios are 
equal due to the high symmetry of the lattice. 
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FIG. 7. The six vanishing ratios of determinants for four different lattices. For comparison 
Ising determinant ratios are shown squared. Coupling constants and aspect ratios for the lattices 
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FIG. 8. The closed differentials Ak,B}~. The differentials are zero on all edges except the ones 
crossed by the respective bk, Ofc loop. These edges are shown in bold on the figure, together with 
the value that the differential takes at that edge. For convenience we use a a>2 loop different but 
homologically equivalent to the one shown in Fig. [lj. 
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FIG. 9. Graphical description of the action of the discrete Hodge operator on the edge values of 
lattice differentials. Four different definitions, related by 7r/2 rotations, are possible. The definition 
given in the text (|5.3| ) corresponds to the drawing A. The four choices have similar properties and 
generate the same period matrix. 
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FIG. 10. The integer sizes characterizing pair of pants decompositions of higher genus lattices 
for genus= 2, 3 and 4. 
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FIG. 11. The Kasteleyn four vertices decoration and the Fisher six vertices decoration for the 
Ising model with four near neighbors. 



34 




FIG. 12. Two dimer configurations: the standard dimer configuration Cq and an arbitrary 
dimer configuration C. The overlap of the two configurations with double edges removed classifies 
C to be of type T(0102) = T(eoee). 
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FIG. 13. An overlap of type T(2001) can be built by overlapping a type T(0001) with a bound- 
ary X and therefore both belong to the same type T(eeeo). 
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FIG. 14. Graphical proof of the existence of 16 independent clockwise odd orientations of the 
genus two lattice. The lattice is redrawn with oriented edge identifications given by the letters 
a to d. The full line gives a spanning tree on the lattice edges were, using the local equivalence 
transformation, we can choose the orientation. On the dashed edges the orientation is determined 
by the clockwise odd condition and on the four rectangular edges the orientation is not determined. 
This gives a total of 2 4 independent clockwise odd orientations. 
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TABLE I. Convergence with lattice size M of the ratios of determinants of the adjacency 
matrices Ri = Pf/Pfe for the critical ISING model in a lattice with w h =0. 66403677026784896368 
w v =0. 20189651799465540849 and aspect ratio (mi,m 2 ,k,ni,n 2 )=(4,2,2,2,4). See section |TTl| for 
discussion. 
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TABLE II. Convergence with lattice size N of the ratios of determinants of the adja- 
cency matrices for the DIMER model in a lattice with z^= 0.780 2^=0.560 and aspect ratio 
(mi,m2,k,m,ri2)=(4,2,2,6,2). Additional discussion is given in section 01. 
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TABLE III. Comparison of ratios of determinants of adjacency matrices with ratios of theta 
functions for the critical ISING model in three different lattices: II, 12 and 13. Each lattice 
is characterized by the couplings and aspect ratio [wh,w v ](mi,m2,k,ni,n2). The L — ► oo ratios 
are compared with theta function ratios for a fitted period matrix O(fit) and the period matrix 
O(eval) evaluated by the procedure of section [fv] (modular choice A). The first two sets of values 
illustrate the fact that re-scaling of the couplings can be compensated by a re-scaling of the overall 
vertical/horizontal aspect ratio. 
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0.13628 
0.170903 
0.171024 
0.000000 
0.828514 
0.035084 
1.000000 



0.136061 
0.000000 
0.000000 
0.000000 
0.000000 
0.828547 
0.000405 
0.692891 
0.000000 
0.136213 
0.170896 
0.171048 
0.000000 
0.828699 
0.03524 
1.000000 



0.136671 
0.000000 
0.000000 
0.000000 
0.000000 
0.827822 
0.000402 
0.691553 
0.000000 
0.136666 
0.171781 
0.171776 
0.000000 
0.827818 
0.035512 
1.000000 



0.218541 
0.000000 
0.000004 
0.000000 
0.000000 
0.815064 
0.000022 
0.596530 
0.000004 
0.133375 
0.270077 
0.184916 
0.000001 
0.729883 
0.051561 
1.000000 



0.218546 
0.000000 
0.000000 
0.000000 
0.000000 
0.815061 
0.000022 
0.596536 
0.000000 
0.133377 
0.270087 
0.184917 
0.000000 
0.729892 
0.051562 
1.000000 



0.219812 
0.000000 
0.000000 
0.000000 
0.000000 
0.814718 
0.000023 
0.594928 
0.000000 
0.133252 
0.271820 
0.185259 
0.000000 
0.728157 
0.052031 
1.000000 



n 



ii 



n 



12 



n 



22 



n 



ii 



n 



12 



n 



22 



n 



ii 



n 



12 



22 



fi(fit) 

ri(eval) 



1.422286 -1.208265 2.416530 
1.42099 -1.20562 2.41124 



1.422067 
1.42056 



-1.205355 2.411027 
-1.20477 2.40953 



1.389893 
1.38931 



-1.297942 
-1.29614 



2.457896 
2.45252 



TABLE IV. Comparison of ratios of determinants of adjacency matrices Ri with ratios of 
theta functions for the DIMER model in three different lattices: PI, P2 and P3. Each lattice is 
characterized by the dimer weights and aspect ratio [zh,z v ](mi,rri2,'k.,ni,n2). 





A: (4,2,2,6,2) 


B: (4,2,2,6,2) with m <-> b { 


C 


!: (2,4,2,2,6) ) 


L 


M 


n u 




^22 






22 


n n 




22 


5 


1000 


1.383663 


-1.279193 


2.401681 


1.414286 


0.7392862 


0.7821439 


1.226958 


-1.122488 


2.401681 


6 


1440 


1.384849 


-1.282726 


2.412281 


1.415022 


0.7407580 


0.7850874 


1.231669 


-1.129555 


2.412281 


7 


1960 


1.385681 


-1.285248 


2.419848 


1.415555 


0.7418236 


0.7872187 


1.235032 


-1.134599 


2.419848 


8 


2560 


1.386311 


-1.287139 


2.425518 


1.415958 


0.7426305 


0.7888325 


1.237552 


-1.138380 


2.425518 


9 


3240 


1.386801 


-1.288608 


2.429926 


1.416274 


0.7432626 


0.7900967 


1.239511 


-1.141318 


2.429926 


10 


4000 


1.387192 


-1.289782 


2.433450 


1.416528 


0.7437711 


0.7911136 


1.241077 


-1.143667 


2.433450 


11 


4840 


1.387513 


-1.290743 


2.436332 


1.416737 


7441889 


0.7919492 


1.242358 


-1.145589 


2.436332 


oo 


OO 


1.38931 


-1.29614 


2.45252 


1.41792 


0.746553 


0.796677 


1.24955 


-1.15638 


2.45252 



TABLE V. Convergence 
h=0.780 v=0.560 and aspect 
homology basis: A, B and C. See section Nl for discussion. 



with lattice size M of the period matrix elements for a lattice with 
ratio (mi,m2,k,ni,n2)=(4,2,2,6,2) for three different choices of the 
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A: (4,2,2,6,2) 


B: (4,2,2,6,2) with a* <-> b { 




C: (2,4,2 


,2,6) 


i 


0;(484O) 


0i(oo) 


i 


0,(4840) 


l (oo) 


i 


(9i(4840) 


0i(oo) 


1 


0.223711 


0.219812 


1 


0.210138 


0.213846 


10 


0.223711 


0.219812 


2 


0.000000 


0.000000 


3 


0.000000 


0.000000 


9 


0.000000 


0.000000 


3 


0.000000 


0.000000 


2 


0.000000 


0.000000 


4 


0.000000 


0.000000 


4 


0.000000 


0.000000 


4 


0.000000 


0.000000 


3 


0.000000 


0.000000 


5 


0.000000 


0.000000 


9 


0.000000 


0.000000 


13 


0.000000 


0.000000 


6 


0.813656 


0.814718 


11 


0.817352 


0.816342 


14 


0.813656 


0.814718 


7 


0.000027 


0.000023 


10 


0.000014 


0.000017 


7 


0.000027 


0.000023 


8 


0.589972 


0.594928 


12 


0.607228 


0.602514 


8 


0.589972 


0.594928 


9 


0.000000 


0.000000 


5 


0.000000 


0.000000 


2 


0.000000 


0.000000 


10 


0.132858 


0.133252 


7 


0.134120 


0.133805 


1 


0.132858 


0.133252 


11 


0.277170 


0.271820 


6 


0.258651 


0.263681 


12 


0.277170 


0.271820 


12 


0.186317 


0.185259 


8 


0.182634 


0.183640 


11 


0.186317 


0.185259 


13 


0.000000 


0.000000 


13 


0.000000 


0.000000 


5 


0.000000 


0.000000 


14 


0.722803 


0.728157 


15 


0.741334 


0.736302 


6 


0.722803 


0.728157 


15 


0.053486 


0.052031 


14 


0.048528 


0.049852 


15 


0.053486 


0.052031 


16 


1.000000 


1.000000 


16 


1.000000 


1.000000 


16 


1.000000 


1.000000 



TABLE VI. The sixteen theta function ratios for the three period matrices calculations of table 
0, both the largest size and the extrapolated value are given. The different ratios are rearranged 
to facilitate comparison. 
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TABLE VII. The relative signs for each dimer configuration type in the various clockwise 



odd edge orientations. The linear combination of the orientations that gives the same sign to all 
configurations is shown in the last row. 
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-Pi 


1 1 

Pi 


^3 


Pa 


P5 


Pq 


P7 


^8 


P9 


P10 




P\2 


P\2, 


P\A 


Pl5 


Pl6 




1 


1 


1 


-1 


1 


1 
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— 1 
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TABLE VIII. The b matrix relating the sixteen partition functions with the Pfafnans of the 
sixteen clockwise odd orientations. Each Zi is given by | of the linear combination of the Pj in 
the corresponding row. 
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